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ABSTRACT 

We investigate the local supermassive black hole (SMBH) density function and relative 
mass accretion rates of all active galactic nuclei (AGNs) identified in a volume-limited 
sample of infrared (IR) bright galaxies (Ltr > 3 x 10 9 L Q ) to D < 15 Mpc (Goulding 
& Alexander 2009). A database of accurate SMBH mass (Mbh) estimates is compiled 
from literature sources using physically motivated AGN modeling techniques (rever- 
beration mapping, maser mapping and gas kinematics) and well-established indirect 
Mbh estimation methods (the M-cr* and MsH-kicbui relations). For the three sources 
without previously published Mbh estimates, we use 2MASS A-band imaging and 
GALFIT to constrain the bulge luminosities, and hence SMBH masses. In general, we 
find the AGNs in the sample host SMBHs which are spread over a wide mass range 
(M B h ~ (0.1-30) x 10 7 M Q ), but with the majority in the poorly studied M B h « 10 6 - 
10 7 Mq region. Using sensitive hard X-ray (2-10 keV) and mid-IR constraints we 
calculate the bolometric luminosities of the AGNs (Lbo1,AGn) and use them to esti- 
mate relative mass accretion rates. We use these data to calculate the volume-average 
SMBH growth rate of galaxies in the local Universe and find that the AGNs hosting 
SMBHs in the mass range Mbh ~ 1O 6 -1O 7 M0 are dominated by optically uniden- 
tified AGNs. These relatively small SMBHs are acquiring a significant proportion of 
their mass in the present-day, and are amongst the most rapidly growing in the local 
Universe (SMBH mass doubling times of w 6 Gyrs). Additionally, we find tentative 
evidence for an increasing volume- weighted AGN fraction with decreasing SMBH mass 
in the Mbh ~ 10 6 -10 8 M range. Overall, we conclude that significant mass accretion 
onto small SMBHs may be missed in even the most sensitive optical surveys due to 
absent or weak optical AGN signatures. 
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1 INTRODUCTION 

It is now well established that all massive galaxies (M» ~ 
10 10 -10 12 Mq) in the local Universe harbour central super- 
massive black holes (SMBHs), with masses proportional 
to those of their stellar spheroids (hereafter, bulg e; e.g., 
iKormendv fc Richstonell 19951 : iMagorrian et al.lll998h . Com- 
parisons between the SMBH mass density in the lo- 
cal Universe and the total energy produced by active 
galactic nuclei (AGNs) across cosmic time have shown 
that these SMBHs were primar i ly grown through mass- 
accre tion events (e.g.. lSoltanll9 82; Rccs 1984; Ma rconi et al.l 
l2004h . The space density of high-luminosity AGNs ap- 
pears to have peaked at higher redshifts than lower- 
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luminosity AGNs, suggesting that the most massive SMBHs 
(M BH » 10 8 -10 9 M ) grew first, a result commonly referred 



to 



as AGN cosmic downsizing' (e.g., Cowie et al.l 20031: 



Ueda et al l 120031; iMcLure fc Dunlod 120041 ; lHasinger et all 



20051 ; lAlonso-Herrero et alj 120081 ) Extrapolation of these 



results imply that the most rapidly growing SMBHs in 
the nearby Universe should be of comparatively low mass 
(Mbh <C 10 8 Mq). To determine the characteristic masses 
of these growing SMBHs requires a complete census of AGN 
activity and SMBH masses in the local Universe. 

Using data f rom the Sloan Digital Sky Survey (SDSS; 
lYork et al.l l2000h in conjunction with the well estab- 
lished SMB H-stellar velocity dispe r sion relation (hereafter , 
M-a*; e.g., iGebhardt et al.l |2000| ; iTremaine etafl |2002| ). 



iHeckman et al. (2004; hereafter, H04) deduced that rela- 



tively low mass SMBHs (Mbh ~ 3 x 10 M ) residing in 
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moderately massive bulge-dominated galaxies host the ma- 
jority of present-day accretion onto SMBHs. However, the 
space density of SMBHs derived from the M— c* relation 
in the optical survey of H04 was limited by the spectral 
resolution of the SDSS (<r„ > 70kms~ 1 ) to SMBHs of 
Mbh i 3 x 10 6 M (assuming the M-cr» relation of Geb- 
hardt et al. 2000). Furthermore, due to attenuation of op- 
tical emission by dust, source selection and AGN classifi- 
cation at optical wavelengths will be biased against gas- 
rich, dust-obscured objects. These surveys are unlikely to in- 
clude galaxies hosting the smallest bulges, and consequently 
the lowest mass SMBHs, and may therefore be missing a 
significant proportion of SMBH growth in the local Uni- 
verse. Indeed, the nearby Scd galaxy, NG C 4945, hosting 
a low-mass SMBH (M BH ~ 1.4 x 10 6 Mm ; iGreenhill et all 
1 19971) only displays evidence for AGN activity in X-ray 
lllwasawa et al.l Il993|) and mid-inf rared (mid-IR) observa- 
tions l|Goulding fc Alexander! l2009to . By contrast, the AGN 
in NGC 4945 (a ccreting at <; 30 percent of the predicted 
Eddington limit; [itoh et al.l [2008^1 is completely hidden at 
optical wavelengths, and classified as a starburst galaxy. 
Clearly, using optical data alone, the intrinsic AGN prop- 
erties of sources similar to NGC 4945 cannot be derived. 



While optical emission-line diagnostics alone cannot re- 
liably characterise the properties of a non-negligable frac- 
tion of the AGN population, they are readily identified at 
obscuration independent wavelengths (e.g., X-ray; mid-IR). 
Hence, the identification of AGNs made at X-ray and mid- 
IR wavelengths complements traditional UV/optical meth- 
ods to yield a more complete census of AGN activity. In- 
deed, using the high resolution mid-IR spectrograph on- 
board the NASA Spit zer Space Telescope (5piteer-IRS), 
iGoulding fc Alexander! (2009; GA09) found using the first 
complete volume-limited sample of all (~ 94 percent) lo- 
cal (D < 15 Mpc) bolometrically luminous galaxies (Lib. > 
3 x 10 9 L©), that « 50 percent of local AGNs are not iden- 
tified in sensitive optical surveys. At least 30 percent of 
these AGNs were previously identified as pure optical star- 
burst galaxies, similar to NGC 4945 (i.e., not even oth- 
erwise known to be tr ansition-type objects as defined by 
iKauffmann et al]|2003h . Furthermore, « 30 percent of the 
optically unidentified AGNs were found to reside in late- type 
spiral galaxies (Sc-Sd; e.g., similar to NGC 4945). Compli- 
mentary to this, from a heter ogeno us sample of Palomar 
galaxies, ISatvapal et al. | |2007h and ISatvapal etafl (|200fet ) 
have also concluded that optically unidentified AGNs exist 
in some late-type spiral galaxies. With the inclusion of these 
new optically unidentified AGNs, it is natural to ask, what 
are the masses of local active SMBHs, what are their Ed- 
dington ratios, and hence, how rapidly are active SMBHs 
growing in the local Universe? 

In this paper, we investigate the growth rates and space 
density of actively accreting SMBHs using the 17 AGNs 
identified in the volume-limited survey of GA09. Whilst 
the source statistics considered here are significantly smaller 
than those studies using the SDSS, this work compliments 
that of H04 by including a relatively large number (given 
the considered small volume) of optically unidentified AGNs 
(10) which would not be reliably identified or characterised 



in the SDSS surveyQ Furthermore, by including a significant 
population of bolometrically luminous (but dust-obscured) 
late-type spiral galaxies (Sc-Sd) we are able to extend the 
SMBH density function to M B h < 3 x 10 6 M . As many of 
the late-type spiral galaxies host small galactic bulges, and 
hence lower mass SMBHs, particular attention is paid to 
obtaining accurate mass estimates for these SMBHs. Given 
their proximity, many of the sources in GA09 are well- 
studied and have multiple estimates of SMBH mass (Mbh) 
from a variety of methods (i.e., reverberation mapping tech- 
niques; mapping of water maser spots; gas kinematical esti- 
mates; the M-cx* relation; correlation of Mbh with the lumi- 
nosity of the galactic bulge); below we discuss the relative 
accuracy of each SMBH mass estimate technique. Further- 
more, to determine the relative mass accretion rates and 
hence average growth times of the SMBHs in our sample we 
require the best available estimates of the AGN bolomet- 
ric luminosity (£bo1,agn)- Here we use two approaches: 1) 
for the AGNs with currently published data, we use high- 
quality well-constrained sensitive hard X-ray (2-10 keV) lu- 
minosities to directly measure LboI.agn; and 2) we accu- 
rately infer Lbo1,agn using a well-constrained hard X-ray to 
high-ionisation mid-IR emission line relation. 

In §2 we outline the construction and basic reduction 
analysis of the AGN sample derived from GA09. In §3 we 
present the SMBH mass estimates. For a minority of ob- 
jects (three out of 17 AGNs) without published Mbh es- 
timates we outline the use of a bulge/disc decomposition 
method with 2MASS if-band images, and following Marconi 
& Hunt (2003), we use the MsH-Lbui relation to estimate 
their SMBH masses. In §4 we use hard (2-10 keV) X-ray 
measurements and high-ionisation mid-IR emission to esti- 
mate the intrinsic luminosity of the AGNs considered in our 
sample. Using our well-defined estimates for SMBH mass 
and AGN bolometric luminosity, we investigate the relative 
mass accretion rates of our sample of active SMBHs in §5. 
We use these estimates to provide new constraints on the 
volume-average SMBH growth rates in the local Universe. 
We furth er compare th e se res ults to the previous works of 
H04 and iGreene fc Hcl (120071 ) by producing a local AGN 
population density function. Finally, in §6 we present our 
conclusions. 



2 THE SAMPLE 

2.1 Selection and Data-reduction 

The sample of local AGNs is derived from the Spitzer-IRS 
spectral investigation of a volume-limited sample of IR- 
bright galaxies to D < 15 Mpc by GA09. High-resolution 
Spitzer-IRS spectroscopy (R ~ 600) was obtained for sixty- 
four of the sixty-eight (~ 94 percent) galaxi es detected in 
the Revised Bright Galaxy Sample (RBGS; ISanders et al.l 



1 Following GA09, throughout this paper we define an optically 
unidentified AGN as an object which is not unambiguously iden- 
tified as a Seyfert galaxy at optical wavelengths using solely tra- 
ditional optical emission-line diagnostics. This has the advantage 
that we may directly compare statistics from the sample consid- 
ered here to those derived from large JV surveys such as the SDSS. 
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with Lir <; 3 x 10 9 L/0 within the considered vol- 
A particular advantage to a volume-limited sample 
is that, unlike magnitude-limited surveys, they do not suf- 
fer from radial selection effects and thus can be used to 
construct volume-averaged statistics. By contrast, volume- 
limited samples yield limited source number statistics due to 
their inability to probe both the faintest and most luminous 
systems. 

Given the relatively small distance scale considered 
here, it is prudent to note the importance of distance 
measurement. For our consider sample, luminosity dis- 
ta nces were ca l culate d using the cosmic attractor model 
of iMould et ail ( 2000) which accurately adjusts heliocentric 
redshifts to the centroid of the local group. Hence, unlike 
many other local surveys, to D < 15 Mpc our sample does 
not include galaxies from local over densities such as the 
Virgo cluster at D ~ 16 Mpc (see Appendix Al for a de- 
tailed analysis and validation of the considered space-volume 
in this survey). 

The Spitzer-IRS spectroscopic data presented in GA09 
was reduced using a custom idl pipeline which utilises 
the Spitzer Science Center data-processing packages SPICE, 
IRSCLEAN and CUBISM. For further detailed information 
on the reduction processes and spectral analyses see §§2.1 
and 2.2 of GA09 and references there-in. 



2.2 [NeV] as an unambiguous AGN indicator 

Due to the very high-ionization potential of [NeV] (97.1 eV) 
we consider its detection coincident with the galactic nucleus 
in mid-IR spectroscopy to be an almo st unambiguous iden- 
tifier of AGN activity. Theoretically, ISchaerer Sz Stasiriskal 
jl999l ) have predicted that extremely hot O and B stars, 
in particular dense populations of Wolf-Rayet stars, may 
produce ionization spectra capable of exciting lines such as 
[OIV] (54.9 eV) and [NeV]. However, observationally, GA09 
found that even in extreme Wolf-Rayet galaxies, whilst 
[OIV] is clearly detected in these types of systems, [NeV] 
emission remains absent to t he detection limit s of this sur- 
vey. Complimentary to this, lHao et all |2009l ) find from a 
Spitzer-IRS study of 12 Blue Compact Dwarf galaxies that 
the mid-IR spectroscopy for eight of their sample contain 
[OIV] emission, however none appear to be producing [NeV]. 

Similarly, extreme starburst driven shocks have also 
been predicted to excite some high-ionization lines 
such as [NeV] (|Allen et al.l l200ct ). however these re- 
quire exceptionally high velocities, and based on the 
[NeV]A14.32 um-[NeII]A12.81 urn and [NeIII]A15.51 um- 
[NeII]A12.81 um emission line ratios presented in GA09, 
the AGNs in this study are not consistent with shock 
models. Furthermore, AGNs which contain strong star- 
formation contributions to their bolometric luminosity are 
often found to have rel atively low mid-I R [NeV] -[Nell] 
and [NeIII]-[NeII] ratios l|Armus et alj|200d : Satyapal et al. 



2 Based on the principle of reprocessed emission in the Unified 
AGN model, selection on the basis of IR-luminosity will select 
all bolometrically luminous AGNs with Lx ^ 10 41 ergs -1 . How- 
ever, we note that it will also select dust-rich star-forming systems 
which may dominate the bolometric luminosity of the galaxy. 



2008; iGoulding fc AlexandeJ [20091 ; iDale et al l [20091 ) . How- 
ever, these ratios are not necessarily strong tracers of the in- 
trinsic power of the AGN. Indeed, many Seyfert 2s are found 
to be hosted in galaxies where-by sta r-formation dominates 
the IR spectral energy distribution (jWeedman et all 120051 ; 
iBuchanan et al.ll2006l ; iDeo et al.ll2007l )7 thus yielding a low 
[NeV]-[NeII] ratio (log-average « 0.02); however, the central 
source may still be extremely luminous at other energies e.g., 
NGC 4945 is the most luminous local AGN at E > 20 keV 
(Done et al. 1996) and by contrast has a [NeV]-[NeII] ratio 
of « 0.01. 

The mid-IR spectra of seventeen of the sixty- four galax- 
ies (« 27^g percent) presented in GA09 were found to con- 
tain the [NeV] A14.32 urn emission line, and hence, host 
AGN activity. These seventeen sources are the main focus 
of the current paper. See Table 1. 



3 BLACK HOLE MASS DETERMINATION 
3.1 Archival Data 

To accurately determine the relative mass accretion rates 
and space density of active SMBHs in the local Universe 
requires reliable SMBH mass (Mbh) estimates for the sev- 
enteen AGNs within our volume-limited sample. Here we 
outline the construction of the heterogeneous database of 
the most reliable available SMBH masses for these AGNs, 
derived from a variety of archival sources and Mbh estima- 
tion methods. 

Many of the sources in the sample are late-type galaxies 
hosting relatively small bulges, and hence low-mass SMBHs. 
In such systems, SMBH mass estimates are often challeng- 
ing to determine as: 1) characteristically low velocity dis- 
persions can be difficult to measure as they are often at the 
resolution limit of published observations, 2) modeling of 
the contamination from composite stellar populations can 
often lead to inconsistencies between published measure- 
ments and 3) SMBH mass relations are poorly constrained 
at Mbh ~ 10 6 M . 

Given the varying degrees of accuracy associated with 
Mbh measurements estimated from differing methodologies, 
we have chosen to prioritise the archival data (which we fur- 
ther expand on here) based upon two broad categories: (i) 
physically motivated AGN modeling techniques (i.e., rever- 
beration mapping, water maser mapping and gas kinematics 
of the central engine) and (ii) indirect estimations from ob- 
servational relations (the M-<r, and MaH-Lbui relations). 
The adopted Mbh for the sources and their associated mea- 
surement methodologies are given in columns 12 and 13 of 
Table 1. The database contains four SMBH measurements 
determined from physically motivated AGN modeling and 
13 from indirect methods. 



3.1.1 Direct SMBH mass constraints from reverberation 
mapping, maser mapping and gas kinematics 

Under the assumption that the gas in the broad-line 
region (BLR) is virialised by the SMBH and the or- 
bital motion of the gas is Keplerian, Mbh estimations 
are possible through rev erberation mapping techniques 
jBlandford fc McKedll982} ). The time-lag between changes 



4 A. Goulding et al. 



ES0121 -G006 













• 


\ 





















-1-2 2 1 -1-2 2 1 
ARA (arcmin) 
06h07m29.5s -61d48m28s 



NGC1448 




1 -1-2 2 1 -1-2 2 1 -1 
ARA (arcmin) 
03h44m31 .9s -44d38m41s 



NGC1792 




05h05m13.7s -37d58m46s 



Figure 1. GALFIT l|Peng et al] I2002T ) two-dimensional 
bulge/disc decompositions for the three AGNs (ES0121- 
G006, NGC 1448 and NGC 1792) with M BH estimated from the 
Mbh- LK,bul relation. The panels show (left) Zf-band Two Micron 
All-Sky Survey (2MASS) extended source image of the galaxy, 
(middle) GALFIT model produced by fitting a Sersic profile, 
an exponential disc and a mean averaged background to the 
observed iC-band image, (right) the residual image resulting from 
the subtraction of the model profile from the observed galaxy 
(i.e., observation — model). See §3.2.1 for a brief explanation of 
the residual images. 



in the AGN continuum flux and the response of the broad- 
line region is used to directly infer the size of the virial 
radius and hence the mass of the SMBH (for an in-depth 
review see Peterson 2001). To date, reverberation mapping 
is widely accepted t o be the most reliab l e of Mbh estimation 
methodologies (e.g . .lWandel et al.lll999l:IPeterson fc Wandell 



a minor- 



1 19991 ; lOnken et alj|2003l ; iBentz et alj 120091 ). Only 
ity of the galaxies in the sample (NGC 1068, NGC 4051 and 
NGC 5033) are known to have detected BLRs. However, the 
BLRs of NGC 1068 and NGC 5033 (both are Seyfert 1.9) 
are extremely weak and are therefore likely to suffer from 
optical extinction. Thus, Mbh estimates using the BLR are 
not possible for these tw o AGNs. Hence, only one object 
(NGC 4051: IWandeill999l ) in our sample has archival rever- 
beration mapping data giving Mbh » 1.41JJ x 10 6 M Q . 

Complimentary to reverberation mapping, and also as- 
suming Keplerian motion, mapping of water maser spots 
and observations of gas kinematics within the gravitational 



sphere of influence of the SMBH are thought to yield rel- 
atively precise measurements of Mbh- Again, due to the 
nature of the observations and the requirement of a suit- 
able gas disc (i.e., a relatively face-on inclination to the ob- 
server), few Mbh estimations using these methods exist in 
the current literature. Indeed, in our sample there are cur- 
rently only two AGNs (NGC 1068 and NGC 4945) with 
Mbh estimations from the mapping of water maser spots 
(Mbh ~ 1.6 x 10 7 and 1.1 x 1O 6 M0, respectively) and one 
AGN (NGC 5128) with a spatially resolved gas dynamical 
Mbh estimate (M BH « 2.4 x 10 8 M Q ). 



3.1.2 Indirect SMBH mass constraints from the M-a* 
relation 

Since the seminal discovery that the mass of the stellar 
spheroid is closely related to Mbh (Magorrian et al. 1998), 
indirect Mbh estimation metho ds have become ubiquitous 
in the current literature (e.g., Ferrarese fc Merrittl | 2000| ; 



Gebhardt et al.ll2000l ; iNelson et al.ll2004l ; lOnken et al.ll2004l ; 
Greene fc Hdl2006l ). Comparisons of spatially resolved stel- 



lar dynamics with reverberation-based Mbh estimates (over 
the range of M BH « 2 x 10 6 -10 9 M ) show that these 
are well-correlated with the effective stellar velocity disper- 
sions (cr*) within the galaxy bulge (i.e., the M-rr* relation; 
iGebhardt et~aT]|200rj| ). 



Mr 



(1.2 ± 0.2) x 10 s M (- 



(3.75±0.30) 



(!) 

.200 km s" 1 J v ; 

with an intrinsic scatter of only 0.34 dex. Furthermore, 
with the inclusion of a sample of dwarf Seyfert 1 galaxies, 
the correlation may be reasonabl y extrapolated to SMBH 
masses in the range 10 5 -10 6 M Q i|Barth et al.ll2005l ). How- 
ever, we acknowledge that to-date, the M-cr* relation is still 
poorly constrained for Mbh < 10 6 Mq and suffers from 
systematic uncertainties. Specifically, these are caused by 
the lack of homogeneous source statistics in this mass re- 
gion and also the poor understanding of the physical nature 
of the central region in very late-type (Scd-Sm) galaxies 
(e.g., the existence of either a pseudo- or classical bulge; 
iKormendv fc Kennicutl |2004| ). Of the three AGNs in our 
sample which are classified as Hubble-type Scd or later, only 
NGC 3621 (Sd) has an Mbh measurement estimated from 
the M-ct* relation0 Specifically, using the Echellette spec- 
trograph on Keck-II, Barth et al. (2009) measure an accu- 
rate line-of-sight stellar velocity dispersion of the Calcium 
II triplet lines (CaT; AA8498, 8542, 8662) observed in the 
central nuclear star cluster. They find a velocity dispersion 
consistent with a SMBH mass of M BH ~ 3 x 10 e M BH . We 
find that all of the AGNs in our sample host SMBHs with 
M B h <i 10 6 Mq and thus are unlikely to suffer significant 
systematic uncertainities arrising from the use of the M-cr* 
relation to estimate Mbh even in the most late-type galaxies. 

Ten of the 13 AGNs in our sample without direct 
Mbh measurements have published Mbh estimates using 

3 The Mbh measurement for NGC 1448 (Scd) is estimated from 
the Mbh-Lk bul relation (see §3.2.1) due to the lack of suffi- 
cient spectroscopic data to measure a stellar velocity dispersion. 
NGC 4945 (Scd) is an H2O megamaser source, and hence has an 
accurate and direct measurement of Mbh (Greenhill et al. 1997). 
See also Table 1. 
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the M-cr, relation. For consistency purposes (and where 
possible) we have used the central stellar velocit y disper- 
sions given in the recently published catalogue of [Ho et al.l 
(2009). They measure the central cr* for t he 486 galaxies 
within the Palomar Survey (|Ho et al.lll997al lbl) using the av- 
eraged val ues derived from th e fitting of stellar absorption 
templates (|Valdes et all l2004h to the blue (4230-51f0 A) 
and red (6210-6860 A) spectral ranges (i.e., the published 
spectroscopy does not include standard velocity dispersion 
measu rement features; e.g., CaT). Where available. iHo et all 
(2009) compare cr* values derived from the modeling of the 
stellar absorption features to previously published measure- 
ments from CaT lines which are available in the HyperLeda 
database Measured errors are compared between the Palo- 
mar cr* measurements and the weighted average adopted 
by HyperLeda for the available published cr* measurements. 
The final adopted measurement of cr, by Ho et al. (2009) 
is that with the smallest overall error. Where the values as- 
sumed by Ho et al. (2009) are previously published or are 
from HyperLeda we quote these references in column 14 of 
Table 1 (6 objects). For the AGNs in our sample which are 
not part of the Palomar survey, values of cr, derived from 
direct fitting analyses of the CaT lines are adopted from 
other published sources (see column 14 of Table 1; 4 ob- 
jects). All final adopted cr, measurements are converted to 
Mbh estimates using Equation 10 

3.2 Galaxy Decompositions using GALFIT 

3.2.1 Indirect SMBH mass constraints from, the 
MBH _ iK,bui relation 

Three of the AGNs within our sample (ESO121-G006, 
NGC 1448 and NGC 1792) currently lack archival direct or 
indirect SMBH mass constraints. Hence, for these three ob- 
jects, we follow the formalism of Marconi & Hunt (2003; 
hereafter, MH03) and use 2MASS _R"-band imaging and 
GALFIT, th e two- dimensional imaging analysis software of 
iPeng et alj l|2002h . to constrain the bulge luminosities and 
therefore, Mbh for these three AGNs. 

Near-IR (0.9-4.8 \an) emission is a strong tracer of stel- 
lar mass and is less susceptible to the effects of dust/gas 
extinction than optical emission. As a result of this, the K- 
band (2.2 um) is shown to provide the strongest correlation 
of all near-IR bands between the luminosity of the bulge and 
Mbh (MH03). 

For the bulge-disc decomposition image analysis, we 
have obtained archival if-band imaging for ESO 121-G006, 
NGC 1448 and NGC 1792. These images were retrieved 

4 The HyperLeda database is a continuously updated elec- 
tronic catalogue of galactic measurements available at 
http: / /leda. univ-fyonl.fr/l Specifically, HyperLeda contains 
a consolidated list of archival velocity dispersions for many 
nearby galaxies. 

5 We note that NGC 6300 curre ntly has two measurements o f 
cr* from fitting of the CaT lines IIGarcia- Rissmann et alj l2005^, 
which were obtained through direct-fitting and cross-correlation 
analyses. For consistency with other measurements, we adopt the 
direct- fitting value of <r* . We find that the derived Mbh from the 
cross-correlation method is a factor two larger; however, using this 
larger Mbh measurement will have little impact on our overall 
results. 



from the Two Micron All-Sky Survey (2MASS) extended 
source catalogue and consist of pre- mosaicked (1 arcsecond 
per pixel resolution) all-sky atlas images. The K-h&nA im- 
ages of the three galaxies were modeled with a central point 
spread function (PSF) and a constant sky background con- 
tribution, whilst the bulge and host-galaxy components were 
modeled using variations of the Sersic profile: 

E(r) = E e e-' ^l(r ' /r ' !)1/ "- 1, (2) 

where r e is the effective radius of the profile, E e is the surface 
brightness at the effective radius, n is the power-law (Sersic) 
index, and k is coupled to n such that half of the total flux 
of the object is within the effective radius. We employ two 
special forms of the Sersic profile in our GALFIT modeling, 
the exponential (n = 1) and the de Vaucouleurs (n = 4) 
profiles, which are classically used to model galactic discs 
and bulges, respect i vely. 

lHaussler et alj (|2007T l have shown that the reliability of 
the fitting parameters produced by GALFIT are strongly de- 
pendent on the initial estimates. GALFIT will, in general, fail 
to find the overall global chi-squared (x 2 ) minimum to the 
fit if the initial estimates are poorly constrained. Thus, to 
reduce this systematic effect, and aid the fitting routine, we 
use a simplified 1-dimensional fit to produce initial estimates 
of the fitting parameters. A 1-dimensional surface-brightness 
slice of the K-b&nd image was taken across the major axis 
of each of the galaxies. A surface brightness profile extend- 
ing from the nucleus was produced by averaging the two 
semi-major axes from the slice, and removing the measured 
background flux. Few spiral galaxies are found to host bulges 
with true de Vaucouleurs profiles, thus a global \ 2 reduction 
process was used to simultaneously fit a generalised Sersic 
profile and a fixed (n = 1) exponential disc to the 1-d surface 
brightness profile. From these, we calculate Sersic and disc 
radii, as well as the Sersic index of the bulge. Combining the 
1-d parameter estimates with the tota l jf-band magnitud e 
from the 2MASS Large Galaxy Atlas (Ijarrett et al.l 120031 ') . 
we generate an appropriate set of constraints and initial pa- 
rameters to be input to GALFIT. 

Using the derived parameter estimates, GALFIT is used 
to fit a generalised 2-d Sersic profile with an exponential 
disc to the Tf-band image. To again aid the GALFIT reduc- 
tion analysis, particular attention is paid to simulating ac- 
curate PSFs for the 2MASS images using known standard 
stars (J. R. Lucey private communication) Results of this 
bulge/disc reduction for the three objects are presented in 
Fig. [T] and column 11 of Table 1. 

We have directly tested our robust GALFIT method us- 
ing the late-type galaxies (i.e., SO or later) presented in the 
dataset of MH03 and find close agreement (« 0.1 dex). We 
do note however, that we find a systematic offset of a factor 
~ 2 in bulge luminosity for the AGNs in MH03 that are 
hosted in low-inclination angle late- type galaxies, which is 
likely to be caused by GALFIT over-estimating the contribu- 
tion of the bulge to the total flux of the galaxy. Indeed, when 
directly comparing a sample of reverberation mapped X-ray 
dete cted AGNs to Mb h estim ations using the MH03 formal- 
ism, IVasudevan et all (|2009l ) find similar results. However, 

B A detailed discussion o f GALFIT problems caused by poor PSF 
modeling can be found in lBentz et all J2009t) 
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Table 1. Catalogue of D < 15 Mpc mid-infrared identified AGNs and derived quantities. 



Common D Hubble AGN log(L [OIV] ) log(L x ) Ref. log(L bol[OIV] ) log(L boljX ) K Tot K Bul log(M BH ) Method Ref. 



Name 


(Mpc) 


Type 


ID 


( erg s 1 ) 


( erg s 1 ) 


b x 


(erg 




( erg s 1 ) 


(mag) (mag) 


(M ) 




M BH 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


E121-G006 


14.5 


Sc 


IR 


39.04 






41 


Si 9 




8.98 


10.91 


6.10+°£ 


LbuI 


14 


NGC 0613 


15.0 


Sbc 


1R 


39.38 


_ 


_ 


42 


26 


_ 


7.03 




7.34+™* 


M-CT* 


15 


NGC 0660 


12.3 


Sa 


IR 


39.71 






42 


69 




7.34 


_ 


7 OC-+0.08 

1 .OO Q -^Q 


M-CT* 


21 


NGC 1068 


13.7 


Sb 


IR,0,X 


41.66 


43.48 


1,9,2 


45 


26 


44.85 


5.79 




7 20+ ' 12 


M 


16 


NGC 1448 


11.5 


Scd 


IR 


39.40 






42 


28 




7.66 


10.64 


5 99+ 11 

J - a -0.52 


LbuI 


14 


NGC 1792 


12.5 


Sbc 


11! 


38.26 






40 


49 




7.01 


9.08 




LbuI 


14 


NGC 3621 


6.6 


Sd 


IR.O 


38.18 






40 


68 




6.60 




6.50±p 3 7 


M-(T» 


22 


NGC 3627 


10.0 


Sb 


IR 


38.38 






40 


95 




5.99 




7 oq+0.10 


M-CT* 


21 


NGC 3628 


10.0 


Sb 


IR 


38.81 






41 


51 




6.07 




D - OJ -0.12 


M-CT* 


20 


NGC 4051 


13.1 


Sbc 


IR.O.X 


39.88 


41.72 


3,2,4,5 


42 


91 


42.71 


7.67 




6 15+ 016 


R 


19 


NGC 4945 


3.9 


Scd 


IR,X 


38.72 


42.49 


6,2,7,4,5 


41 


40 


43.61 


5.23 




6 04+ - 05 


M 


17 


NGC 5033 


13.8 


Sc 


IR,0,X 


39.08 


40.85 


8,9 


41 


86 


41.73 


6.96 




7 62+ - 09 
'• DZ -0.16 


M-cr, 


21 


NGC 5128 


4.0 


SO 


IR,X,R a 


39.38 


41.85 


2,9,4,5 


42 


26 


42.86 


3.94 




fi qo+0.20 
8 - da -0.2fi 


G 


18 


NGC 5194 


8.6 


Sbc 


IR,0,X 


38.85 


41.11 


10,2 


41 


56 


42.00 


5.92 




6.881°;" 


M-cr„ 


15 


NGC 5195 


8.3 


Irr 


IR 


37.89 






40 


30 




6.25 




7 ni+0.07 
' • J -0.13 


M-CT* 


20 


NGC 5643 


13.9 


Sc 


IR.O.X 


40.43 


41.08 


11,12,2 


43 


63 


41.98 


7.17 




6.44l°-ii 


M-cr„ 


23 


NGC 6300 


13.1 


Sb 


IR,Q,X 


39.78 


41.63 


13,4,5 


42 


79 


42.60 


6.93 




6.80t°£ 


M-cr, 


21 



notbs: (1) Common galaxy name. (2) Distance to source in megaparsecs f rom the Revised Bright Galaxy Survey (RBGS; Sanders 
et al. 2003). (3) Morphological type from RC3 (jde Vaucouleurs et ail 1 199 li b (4) Waveband of AGN identification; IR: Mid-Infrared 
spectroscopy (GA09); O: Optical spectroscopy (references presented in GA09); X: X-ray spectroscopy (2-10 keV; see column 7 for 
references); R: Radio observations. (5) Logarithm of [OIV] A25.89 um luminosity in ergs - calculated using [OIV] flux presented in 
GA09: mean uncertainty is approximately 10 percent. (6) Logarithm of absorption corrected hard X-ray luminosity (2-10 keV) in ergs -1 
which have been converted to the distances given in column 2. (7) Reference for X-ray data. (8) Logarithm of bolometric luminosity 
of the AGN estimated from £[oiv] using Equation 4. (9) Logarithm of bolometric luminosity of the AGN estimated from Lx,2— lOkeV 
using the bolometric corrections described in Marconi et al. (2004). (10) Total J\-Band magnitu de from 2MASS Large Galaxy Atlas 
(Jarrett et al. 2003). (11) i\-Band magnitude of bulge produced using GALFIT llPeng et al.ll2002l ; see section 3.2.1). (12) Logarithm of 
estimated black hole mass and associated 1-sigma errors in solar masses. (13) Method of Mbh measurement; M: Maser Mapping; G: 
Gas Kinematics; R: Reverberation Mapping; M-cr*: Mass- Velocity Dispersion Correlation; Lb u i: If -band Luminosity-Bulge Correlation. 
(14) Reference for Mbh measurement. 

(1 1 iMatt et alj <1997t); (2) iDadinal feOOTf); (31 iPounds et all [I2004T); (41 iTueller et all ll2008h; (51 IWinter et all j2009h : 



REFERENCES: 



(61 iGuainazzi et alj feOOd); (71 lltoh et al, 
iMaiolino et alj dl998ft; (12llG uainazzi et al 



1200811 ; (81 ICappi et al.l <200d); (9 1 iBird et al.l [120071) ; ( 101 iFukazawa et alj feOOll); (111 
120041) ; (13l lMatsumoto et aljfeoMI; ( 141 T his Paper; (15) H vper Leda; (16llGreenhill et al] 



1996h; (17l lGreenhill et al.llll997ll; (18 1 iMarconi et alJfeOO ll); (191 IWandell jl99sh : (201 iHo et all j2009h; (2 1l lBarth et alj [120021) ; (221 



iBarth et ail l l2009t) ; (23) IWhittlj i ll992l) ; (241 iGarcia-Rissmann et alj d2005h ; a for a review see llsraell lll998ll 



the three galaxies fitted in our sample are all moderately to 
highly-inclined and thus this systematic effect will be negli- 
gible. 

In Fig. [T] we show the three GALFIT produced image 
cubes obtained following our bulge/disc fitting routines. 
Within each of the residual (observed - model) images it 
is clear that the bulge is well fitted by a Sersic profile. 
The edge-on galaxy ESO121-G006 is well fit by an expo- 
nential disc combined with a Sersic profile, with no dis- 
tinguishing residual features. The residual of NGC 1448 
highlights the existence of its spiral arms and shows the 
presence of a truncated disc combined with a possible bar 
structure which our simplified modeling technique is inca- 
pable of fitting; however, the bulge fit does not appear to 
be compromised. Indeed, our derived Mbh estimations for 
ESOT21-G006 and NGC 1448 are consistent with the Mbh 
uppe r limits obtained from ste llar mass-to-light ratio anal- 
yses (|Ratnam fc Saluccill2000l ). The residual image of the 
moderately inclined galaxy, NGC 1792, contains strong spi- 
ral arms as well as a point-like nuclear source. We note that 



due to the inclination angle of this source, the derived Mbh 
may be systematically over-estimated by a factor of « 2 (see 
above). Thus, our derived Eddington rate for NGC 1792, 
presented in §5.1, should be considered a lower limit. 

From the obtained A"-band bulge magnitudes (column 
11, Table 1) we calculate bulge luminosities (Z/K.bul). Us- 
ing the well established Mbh to bulge _ft"-band luminosity 
relation of MH03 (hereafter, MBH - £K,bui), 

log Mbh = (8.08 ± 0.10) + (1.21 ± 0.13)log ( 8 ^ff ) (3) 

we calculate Mbh for ESO 121-G006, NGC 1448 and 
NGC 1792. The 1-a uncertainties associated with the Mbh 
values are assessed by combining in quadrature the bulge 
magnitude errors as calculated by GALFIT with the intrinsic 
dispersion observed in the MBH - LK,t>ui relation. The resul- 
tant Mbh estimates and 1-a uncertainties are provided in 
column 12 of Table 1. 
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4 BOLOMETRIC CORRECTIONS 



To determine the relative mass accretion rates, and hence 
the average growth times of SMBHs in the local Universe we 
require relatively accurate estimates of the bolometric lumi- 
nosities of the AGNs ( £bo1,agn) in our sample. We use two 
methods to estimate LboI.agn: 1) a direct approach using 
the best available measured hard X-ray (2-10 keV) lumi- 
nosities, and 2) a well-constrained I/bo1,agn _ [OIV] luminos- 
ity relation to infer the intrinsic lum inosity of the AGN (e.g., 
iDasvra et afll200Sl : iMelendez et all 2008. hereafter, M08). 



4.1 Hard X-ray Luminosity as a Tracer of the 
Bolometric Luminosity of an AGN 

High quality hard X-ray spectral analyses arguably provide 
the most unambiguous method for measuring the intrinsic 
luminosity of an AGN since: 1) X-rays are relatively un- 
affected by dust extinction; 2) intrinsic absorption can be 
directly constrained from high S/N data; 3) star-formation 
contamination is often found to be negligible. We therefore 
divide our sample into two categories based on the quality 
and energy range of their available published X-ray data: 
1) AGNs with high S/N X-ray spectra where Nh has been 
accurately constrained and/or AGNs with E > 10 keV 
constraints where the observed X-ray emission will only 
be strongly absorbed for heavily Compton-thick (Ah > 
10 25 cm" 2 ) sources (8 AGNs); and 2) those AGNs with no or 
low S/N X-ray data, i.e., where there are insufficient counts 
to accurately determine Ah and the X-ray flux could have 
large contributions from star- formation (9 AGNs). Hence, 
we specifically do not estimate LboI.agn for those AGNs 
with Lx measurements using Chandra that also have no 
further hard X-ray spectral constraints (E > 10 keV) due 
to the limited band-pass of the instrument at z ~ (0.5- 
8 keV). 

For the AGNs in our sample currently with either high 
S/N spectroscopy or E > 10 keV constraints, we estimate 
£bo1,agn using the Eddington ratio independ ent AGN bolo- 
metric corrections ou tlined in Equation 21 of Marconi et al.l 
(12004 ). We note that I Vasudevan fc Fabianl l|2009r ) have sug- 
gested that the bolometric correction factor (K2-iokcV = 
iboi/^2-iokov) may be a function of the Eddington ratio of 
a considered source. Values of K2-iokcv ~ 10-30 are con- 
sidered to be relatively low bolometric corrections and are 
generally fou nd in AGNs with n < 0.1 (e.g., Vasudevan & 
Fabian 2009; IVasudevan et alj|2010h . For the sample of X- 
ray detected AGNs considered here, we calculate similarly 
consistent values of K2-iok c v ~ 8-30, and thus conclude that 
the Eddington ratio is unlikely to be dominating the bolo - 
metric corrections adopted here from lMarconi et all |2004T l. 

For further consistency, all archival X-ray luminosities 
were adjusted to the distances adopted in column 2 of Ta- 
ble 1. Final adopted Lx measurements, estimated LboI.agn 
from Lx , and archival references for Lx are given in columns 
6, 7 and 9 of Table 1, respectively. 



4.2 [OIV] Luminosity as a Tracer of the 
Bolometric Luminosity of an AGN 

For those galaxies in the sample without good-quality hard 
X-ray constraints we require an additional approach to es- 
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Figure 2. X-ray derived AGN bolometric luminosity 
(^Boi.AGN.x) versus [OIV] A25.89 urn luminosity (£[orv])- 
The objects shown are from this work (blue squares) and those in 
the Swift-HAT survey with detected [OIV] emission (M08; grey 
squares). The solid line is a best-fit bisector and was obtained 
using the iDL-based robustxinefit algorithm. The dashed 
line is the £bo1,AGN~^[OIV] relation of Dasyra et al. (2008) 
derived from the optical luminosity at 5100A. We additionally 
highlight those galaxies in M08 with inconsistent 2-10 and 
14-195 keV luminosities (open circles). These seven galaxies and 
their alternatively published X-ray luminosities (which we adopt 
in our analysis; dotted line) are shown in Table 2; we find that 
five are identified as Compton thick AGNs in the literature and 
the other two are X-ray variable. The galaxies from our D < 15 
Mpc sample which lie significantly off the derived correlation are 
labeled (NGC 4945 and NGC 5643); see text for further details. 



timate the intrinsic luminosity of the AGN, and hence 
the relative mass-accretion rate. Here we build-upon an 
LboI,agn estimation which relies on the AGN-produced 
[OIV] A25.89 urn luminosities (Lrotvi) f° r the AGNs in our 
sample (e.g., Dasyra et al. 2008; M08). 

Based on the simplest Unified Model of AGN 
l|Antonuccilll993| j. the hot dust within the predicted torus, 
close to the central engine, reprocesses absorbed UV, optical 
and X-ray emission into mid-IR emission. Hence, AGN emis- 
sion detected at IR wavelengths is likely to be isotropic and 
independent of viewing angle. As discussed in §2, the detec- 
tion of high-ionisation [NeV] emission (97.1 eV) coincident 
with the nucleus of a galaxy is considered a robust indica- 
tor of AGN activity (e.g., Armus et al. 2006; GA09). Com- 
plimentary to this, GA09 find that [NeV] emission is also 
well correlated with [OIV] emission (54.9 eV) with an in- 
trinsic scatter of only 0.24 dex. As [NeV] emission, and thus 
[OIV] emission, do not suffer from significant star-formation 
contamination and are both comparatively extinction-free, 
they may be used as relatively clean proxies for the bolo- 
metric luminosity of the AGNs presented here (£bo1,agn)0 



7 We note that since star-formation can also produce [OIV] emis- 
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Table 2. Catalogue of revised 2-10 keV luminosities for a subset 
of M08 sample. 



Common 


log(L 2 _i0kcV,M08) 


log(L 2 _lOkcV,alt. ) 


C-thick 


Ref. 


Name 


(ergs" 1 ) 


(ergs" 1 ) 


AGN? 




(1) 


(2) 


(3) 


(4) 


(5) 


Circinus 


40.58 


42.04 


V 


1 


Mrk 3 


41.95 


43.20 


V 


2 


NGC 1365 


40.99 


42.40 


V 


3 


NGC 2992 


41.69 


42.50 


X 


1 


NGC 3079 


40.02 


42.18 


V 


5 


NGC 4388 


41.91 


42.57 


X 


6 


NGC 6240 


42.23 


44.00 


V 


7 



notes: (1) Common galaxy name. (2) Logarithm of 2-10 keV lu- 
minosity adopted by M08. (3) Logarithm of absorption corrected 
2-10 keV luminosity adopted in this work from individual stud- 
ies (2-10 keV luminosities were adjusted using our adopted dis- 
tances). (4) Is the AGN a Compton thick source? (5) Reference 

for adopted 1/2— inkeV measu r emen t. 

references: (11 lYane et al.l J2009h; (2)lAwaki et alj 11200811; f3) 
iRisaliti etaLI j|2009l1; (4 l lYaqpob et al.l 1120071); (5 1 llvomoto et al.l 



j200lll ; (6l lShirai et al.l d2008h ; f7l lVignati et all l|l999ll 

Indeed, for a sample of 35 well-studied optically unobscured 
AGNs, Dasyra et al. (2008) show that both [NeV] and [OIV] 
emission are well correlated with the luminosity of the 5100 
A optical continuum, and hence Lbo1,agn with an intrinsic 
scatter of 0.46 and 0.47 dex, respectively. However, the re- 
lation of Dasyra et al. (2008) is derived from AGNs with 
L [0 iv] > 2 x 10 40 ergs' 1 (XpveV] > 7 x 10 39 ergs' 1 ), and 
hence we test whether it may be reliably extrapolated to the 
more modest luminosity AGNs considered here (log-average 
L [0 iv] ~ 2 x lO^ergs" 1 ). 

We combine our robustly adopted I/boI.agn from 2-10 
keV flux measurements with the X-ray catalogue of nearby 
(z < 0.08) Seyfert galaxies in the Swift-HAT survey which 
have published [OIV] luminosities in M08. The catalogue of 
sources in M08 contain 2-10 keV luminosities obtained pri- 
marily from ASCA data, 14-195 keV luminosities from the 
Swift-TSAT survey and [OIV] luminosities from Spitzer-IRS 
spectroscopy. Combining the M08 sample with our 8 AGNs 
with high-quality X-ray constraints, the range covered in 



L 



[OIV] 



(0.7-7000) x 10 39 ergs 1 . For consistency, we 



convert the M08 2-10 keV luminosities to I/boI.agn using 
the same bolometric corrections adopted in §4.1. In Fig. [2] 
we plot LboI.agn versus L[oiv] f° r the M08 sample and the 
8 AGNs in our GA09 sample with good-quality hard X-ray 
constraints (grey filled squares and blue filled squares, re- 
spectively) . 

Seven of the AGNs in M08 were found to have unusually 
low 2-10 keV luminosities when compared to the quoted 14- 
195 keV luminosities, which are highlighted with open circles 
in Fig.[2]and are shown in Table 2. Using high-quality X-ray 
spectral analyses, in the literature we find that five of these 
AGNs are identified as Compton thick AGNs (see column 4 



sion, the [NeV]— [OIV] relation may be unreliable for sources with 
exceedingly high star- formation rates (e.g., Ultra Luminous In- 
frared Galaxies with Ljp > 10 12 ergs -1 ; ULIRGs). However, 
there are no ULIRGs within the sample considered here. 



of Table 2). This additional obscuration does not appear to 
have been accounted for in the adopted value of M08, and 
therefore we have chosen to select measurements of Z^-iukcV 
from the literature, which now give good agreement with 
14-195 keV luminosities (see columns 2, 3 and 5 of Table 2; 
dotted- lines in Fig.H}. The other two AGNs (NGC 2992 and 
4388) in M08 with inconsistent 2 -10 and 14-195 keV fl uxes 
are found to be highly variable l|Beckmann et alll2007i and 
lElvis et alj|2004 respectively). 

With the inclusion of the obscuration-corrected 2- 
10 keV luminosities to infer LboI.agn, we find a strong cor- 
relation between LroiV] an d ^BoI.agn which is characterised 
by the equation: 



lop 



10 44 ergs" 1 
(1.31 ± 0.09)log 



(0.38 ±0.09) + 



-^[OIV] 

10 41 erg s 



(4) 



with an intrinsic scatter in the data of ~ 0.35 dex. We 
find good agreement with the Dasyra et al. (2008) relation 
(dashed-line) in the region L[oiv] ~ 10 40 -10 4 



erg s 



where the Dasyra et al. (2008) relation is well sampled). 
However, for -L[oiv] < 10 40 ergs -1 we show that the Dasyra 
et al. (2008) relation will over estimate Lbo1,agn by typically 
« 0.5 dex. 

In Fig. [2] we highlight two AGNs from the D < 15 Mpc 
sample (NGC 4945 and 5643) which appear to be significant 
outliers of the observed correlation. NGC 5643 possibly har- 
bours a variable central source. F r om th e detection of strong 
Fe Kq emission, iMaiolino et al.l (|l998l ) suggest from using 
BeppoSAX data that NGC 5643 is possibly Compton thick 
(Nh > 10 25 cm -2 ); however, Guainazzi et al. (2004) find us- 
ing XMM-Newton data, that it may be Compton thin with 
Nh ~ (6-10) x 10 23 cm -2 . Thus, from current available data, 
the true intrinsic luminosity of the AGN is highly uncertain. 
Here we conserv a tively adopt the Compton thin Lx value of 
iGuainazzi et all ll2004t) ; howe ver, we note that if we use the 
value of IMaiolino et all l|l99^ , then NGC 5643 would lie on 
our derived relationship. NGC 4945, by contrast, has obser- 
vations using Ginga, ASCA, OSSE, the Rossi X-ray Timing 
Explorer, BeppoSAX , Suzaku and most recently Swift-TSAT 
lllwasawa et al.lll993l . iTanaka et all 1 1994 iDone et al.lll996l. 



2008 



iMadeiski et al.ll2000l. IGuainazzi et al.ll2000l . Iltoh et al 
and lTueller et aJj|2008| , respectively) all of which provide ex- 
cellent and consistent spectral constraints over a wide X-ray 
band (1-200 keV) showing that the central source is Comp- 
ton thick (N H ~ 4 x 10 24 cm~ 2 ), with Lx,2-iok c v ~ (3- 
10) x 10 42 ergs _1 and an observed intrinsic variability of a 
factor m 2. Here we adopt the luminosity from the most 
recent observation by Swift, Lx ~ 3 x 10 42 ergs _1 . On 
the basis of our mid-IR constraints, it would therefore ap- 
pear that NGC 4945 is under-luminous in [OIV] flux by 
a factor of w 50 (see Table 1 and Fig. [5J. We suggest 
this deficit in observed [OIV] flux is unlikely to be due 
to host galaxy extinction; the required absorption to ac- 
count for a factor of ~ 50 flux difference is Ay ~ 240 mags 
(Nh ~ 5 x 10 23 cm -2 assuming typical dust-to-gas ratios; 
using Ay /E(B — V) = 3.1). Another possible explanation is 
a temporary decoupling of the X-ray emitting and narrow- 
line regions (i.e., the highly-luminous state of NGC 4945 may 
be a somewhat recent event). Given the spatial difference 
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Figure 3. AGN bolometric luminosity (I/Bol AGNi m ergs -1 ) is plotted against SMBH mass (Mbh) f° r the D < 15 Mpc mid-infrared 
identified AGNs presented in GA09. Associated 1— a error bars for Mbh and Lbo1,AGN estimations are shown (see §3 and §5.1, respectively 
for details of their derivations). AGNs which are previously identified in optical surveys are highlighted with open circles. AGNs with Mbh 
estimates from reverberation mapping (downward triangles), gas dynamics (hour glass), maser mapping (diamond), the M— cr* relation 
(squares) and the Mbh-Lk bul relation (upward triangles) are plotted. NGC 5195 is represented with both an upward triangle and 
square (with a dashed-line connector) as both of the Mbh estimates for this galaxy are highly uncertain given its irregular morphology. 
Constant ratios of Eddington luminosity and their implied SMBH mass-doubling times are illustrated for ij = 10~ 3 , 1CT 1 , 1 (30 Gyrs, 
300 Myrs and 30 Myrs; solid line, short-dash and long-dash, respectively). Contours are shown for the active galaxies in the SDSS optical 
survey of H04. In general, we probe lower SMBH masses and AGN luminosities than those of H04, and we find that the majority of these 
AGNs would not be detected using optical SDSS data alone. 



(and hence, the light-travel time) between the two emission 
regions: 10~ 4 pc and 1-10 pc, respectively, the photoionisa- 
tion of the narrow-line region, and thus the observed [OIV] 
emission, may take w 100 years to respond to the changes in 
the X-ray emitting region. We therefore suggest that the in- 
trinsic scatter in the observed [OIV]-Lboi,AGN relation may 
be significantly reduced if it was possible to account for vari- 
ability in the central region of all of the AGNs. Indeed, we 



find the average dispersion decreases to w 0.31 dex if we 
remove NGC 4945 from our analysis. 

Using Equation 4 and L[oiv] from GA09 (column 5 of 
Table 1), we estimate -£/bo1,agn 

(column 7 of Table 1) for 
those AGNs in our sample currently without good hard X- 
ray measurements. We use these LboI.agn estimates to as- 
sess the relative mass accretion rates (r? ~ Lbo1,agn /LEdd) 
of the SMBHs in our D < 15 Mpc sample. 
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5 RESULTS AND DISCUSSION 

From a volume-limited sample of 64 bolometrically lumi- 
nous galaxies to D < 15 Mpc («j 94 percent complete) 
GA09 unambiguously identified seventeen (« 27^g per cent) 
sources to be hosting AGN activity in galaxies with Lir > 
3x 10 9 L , using [NeV] A14.32 vim emission as a robust AGN 
indicator. Using the SMBH mass and AGN bolometric lumi- 
nosity estimates derived here, we discuss the relative mass 
accretion rates of these seventeen AGNs and use them to de- 
rive the average present-day growth times of SMBHs in the 
very nearby Universe. Furthermore, we evaluate the unique 
contribution that our new optically unidentified AGNs make 
to the space density of active SMBHs in the local Universe, 
which have until now been previously derived from large- 
scale optical surveys (e.g., H04; Greene & Ho 2007). 

5.1 Derived AGN Properties and Relative Mass 
Accretion Rates 

In Fig. [3] we plot LboI.agn against our adopted Mbh es- 
timates (the associated 1-a errors for Mbh measurements 
are described in section 3) for the 17 AGNs in our volume- 
limited sample. LboI.agn is inferred from either accurate 
intrinsic high-quality hard X-ray (2-10 keV) constraints 
(where available) or AGN-produced [OIV] A25.80 urn emis- 
sion (see §4.1 and 4.2). The Z/boI.agn 1-c errors for the 
sources with hard X-ray constraints are the result of com- 
bining the uncertainty in the Lx,2~iokeV measurement with 
that of the mean s pread in the bolome tric correction fac- 
tor employed from iMarconi et all (|2004r ). For those AGNs 
with Z/Bo1,agn derived from I/[oiv]i the error is derived from 
the uncertainty in I/[oiv] as quoted in GA09 combined in 
quadrature with the intrinsic scatter of the empirical [OIV]- 
£t>oi,AGN relation (equation 4). 

Seven objects within the sample have been classi- 
fied as optical AGNs from previous surveys (see GA09 
and Table 1) using typical optical emission-line diagnos- 
tics (e. g., the Baldw in-Phillips- Terlevich diagnostic dia- 
grams; iBaldwin et al.l 1981); however, all have detected 
[NeV] A14.32 um emission, and thus are unambiguously 
identified to host AGNs at mid-IR wavelengths (GA09)Q 
We find that with the exception of NGC 5128 (Centau- 
rus A), our sample is dominated by AGNs with SMBHs 
in the mass ranee M BH ~ (0.1-5) x 10 7 M (median of 
Mbh ~ 7 x 10 6 Mq). Due to the irregular structure of one of 
the galaxies in the sample (NGC 5195), Mbh is poorly deter- 
mined; in Fig. Owe plot Mbh estimates from both the M-<t„ 
and MflH-LK.bui relations (connected blue-dashed line). 

We find the AGNs in our sample are spread over a 
wide range of bolometric luminosities, LboI.agn ~ 10 40 - 
10 45 ergs _1 . To assess the relative mass-accretion rates of 
the sample (I/Boi,AGN/£Edd ~ rj), we over- plot lines of con- 
sta nt Eddingto n ratios (rj ~ 10 -3 , 10 , 1.0; derived follow- 
ing iRees! 1 19841 ) and their associated mass-doubling times 
(t ~ 30,0.3,0.03 Gyrs, respectively). Given the large range 

8 We note that NGC 3627 is ambiguous ly classified as T 2/S2 
from the optical spectroscopy presented in lHo et al.l dl997al') . and 
thus m ay host an active central source. However, I Roberts et al.l 
112001! ) suggest from its optical emission-line ratios that NGC 3627 
is most likely a LINER/HII composite. 
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Figure 4. Mbh is plotted against the characteristic mean mass 
doubling time (<2m) of a SMBH in units of Hubble-time (tn) 
for the D < 15 Mpc AGNs. Growth time errors are calculated 
from the log-normal standard deviations of the sample. For com- 
parison, the growth time function of H04 is shown (dashed line). 
We find good agreement with H04 over the region Mbh ~ 0.3— 
10 X 10 7 Mq and further extend the growth time constraints to 
lower SMBH masses (M BH ~ 10 6 M ). 



in bolometric luminosities, it is not surprising that the AGNs 
in the sample are found to be accreting at rates covering over 
5 orders of magnitude (rj « 10 _5 -1). With the exception of a 
few AGNs, the observed range in Eddington ratios is found 
to be roughly consistent with those found by H04 for active 
galaxies (solid contours in Fig. (3J . 

As our work is not limited by the spectral resolution 
of the SDSS (i.e., with a limit of Mbh > 3 x 10 6 M ), we 
show in Fig. Othat significant accretion, rj > 10~ 3 (i.e., ra- 
diatively efficient accretion systems; e.g., thin discs) occurs 
onto SMBHs with Mbh ~ (1-3) x 10 6 M Q . The majority of 
these low-mass, rapidly-accreting SMBHs are hosted in late- 
type, disc-dominated spiral galaxies (Sc-Sd). By contrast, it 
is generally assumed that gas-rich late-type spirals are pref- 
erentially inactive galaxies and that a large bulge may be a 
necessary component for the existence of a SMBH, and thus 
a luminous AGN. Furthermore, of the four AGNs within the 
sample with SMBHs consistent with Mbh ~ 10 6 Mq , we find 
that three sources are not identified as AGNs in sensitive op- 
tical surveys. This indicates that significant SMBH accretion 
may be missed by statistically-large optical surveys such as 
H04 even if the spectral resolution was sufficient to identify 
SMBHs down to Mbh ~ 10 6 Mq. 

For the subset of our AGN sample which host SMBHs 
with Mbh ^3x 10 6 M , we find that many of the optically 
unidentified AGNs are accreting at relatively low Eddington 
ratios (rj iS 10 ), and are unlikely to make a significant ad- 
ditional contribution to the present-day growth of SMBHs. 
However, these same AGNs may form part of a separate, un- 
derlying population of radiatively inefficient accretion sys- 
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tern s such as advection d ominated accretion flows (ADAFs; 
e.g., iNaravan fc Yil 1 19941 ) or those which contain optically- 
thick slim-discs. Further spectral analysis of the X-ray data 
may distinguish between these particular accretion systems, 
but is beyond the scope of these analyses (see Goulding et al. 
in preparation). 

5.2 The Present-Day Growth of SMBHs 

Using the relative mass accretion rates estimated for our 
sample (Fig. [3]), we can infer the volume-averaged growth 
time of SMBHs in the local Universe. Assuming a mean 
Kerr spin parameter (a) for our sample of a ~ 0.67 (e.g., 
iTreister fc Urry||2006l ; lHopkins et alj2007t ). i.e., an accretion 
efficiency (e) of w 0.1, the characteristic mass doubling time 
(fajkr) of a SMBH accreting matter at the Eddington limit is 
tiM ~ 30 Myrs (Rees 1984). Under the further assumption 
that o, and hence e , does not vary significantly for changes 
in Mbh (|King et al.| [200ct) , we assess the present-day growth 
rate of SMBHs0 Following H04, we calculate and extend to 
lower masses (Mbh < 3 x 10 6 Mq) the integrated growth 
of SMBHs. Growth time errors are calculated from the log- 
normal standard deviations of the sample. We note here that 
we also include the optically unidentified AGNs which would 
not be detected in the SDSS. 

In Fig. [4l we find that the mean growth time for low- 
mass SMBHs (M BH « 1O 6 M ) is « 6±g Gyrs, which is 
consistent with these AGNs growing on time-scales similar 
to that of the age of the Universe. Our results are found 
to be broadly consistent with a simple extrapolation of the 
growth times calculated by H04 to Mbh ~ 10 6 Mq (dashed- 
line in Fig. [4j. Thus, the AGNs hosting SMBHs in the mass 
range Mbh ~ 10 6 -10 7 Mq , which are dominated by optically 
unidentified AGNs (see Fig. [3}, are acquiring a significant 
proportion of their mass in the present-day, and are amongst 
the most rapidly growing SMBHs in the local Universe. Fur- 
thermore, we find our derived growth times of SMBHs with 
Mbh <; 3 x 10 6 Mq are in good agreement with those pre- 
sented in H04, with mean growth times of <2M ~ 47^g 
and « 198tgf Gyrs for AGNs with Mbh ~ (0.3-3) and 
(3-30) x 1O 7 M , respectively^ 

5.3 Space-Density of AGNs in the local Universe 

An accurate active SMBH mass function, especially for lower 
mass SMBHs (Mbh ~ 10 6 Mq), is crucial for extending our 
understanding of the role played by accretion in the growth 
of SMBHs across cosmic time. In this section we calculate 
the space density of active SMBHs for our sample and com- 
pare it to complimentary optical studies of local narrow-line 
(NL; H04) and broad-line (BL; Greene & Ho 2007) AGNs, 

9 We note that the spin variation and spin directionality of 
SMBHs in AGNs is currently an ongoing area of research, 
and a consensus between groups has yet to be reached for 
an averag e value of the Ker r spin parameter; for e xam- 
ple see iBrenneman fc Reynolds! J2006I) . iKing et al] i2008h and 
iFabian et al.1 ^2009^ . 

10 We note that our sample contains only one galaxy with Mbh >§; 
10 s Mq (NGC 5128) and thus may not be representative for high 
Mbh systems. 



and the total mass function of local SMBHs by Marconi 
et al. (2004). 

Following Greene & Ho (2007), in the top panel of Fig. [5] 
we plot the volume-weighted space density, $ against Mbh 
in mass bins of 0.5 dex. The volume, V, encompassed by 
the GA09 sample to D < 15 Mpc is V « 1.3 x 10 4 Mpc 3 (3 
As we note in §2, given the luminosity limit imposed in our 
volume-limited survey (Lir ii 3 x 10 9 Lq), we do not include 
the dwarf like systems which are likely to host the very small- 
est SMBHs (Mbh < 10 5 M Q ). With an adjustment for the 
distance model a dopted in this w ork, an examination of the 
Palomar survey (|Ho et al.lll997al ) shows there are possibly 
three optical Seyferts (NGC 185; NGC 1058; NGC 4395) 
with Mbh < 10 J Mq which are not included in our sample 
due to our lower luminosity limit. However, we also note 
that NGC 185 has sin ce been re-classified as an HII galaxy 
|Ho fc Ulvestadll200ll ). 

The derived volume-weighted space density for our ac- 
tive SMBHs (filled squares), which is dominated by NL- 
AGNs, is found to be significantly greater (a factor of 
« 100) than the SMBH density of BL-AGNs (filled cir- 
cles) presented in Greene & Ho (2007) in the mass region 
Mbh ~ (0.9-90) x 10 6 Mq. The significant increase in active 
SMBH density when compared to the BL-AGN density of 
Greene & Ho (2007) is partially to be expected due to the 
greater relative sensitivity of our Spitzer-IRS observations 
coupled with the greater abundance of observed Seyfert 2 
to Seyfert 1 galaxies identified in the local Universe. How- 
ever, this still may not b e a good indicato r of th e intrinsic 
Seyfert l:Seyfert 2 ratio. iTommasin et al] |2010T l find in a 
large sample of 81 Seyfert galaxies that only six sources do 
not contain significant [NeV] A14.32 [i.m emission in their 
mid-IR spectroscopy, the majority of which are Seyfert Is. 
It is likely that the identification of low equivalent-width 
emission lines (such as [NeV] or [OIV]) in BL-AGNs is fur- 
ther complicated by a strong IR continuum emission which 
dominates the mid-IR regime. Hence, it is possible that by 
requiring the detection of [NeV] A14.32 um to infer AGN 
status, we may be rejecting broad-line objects, and thus find- 
ing a lower Seyfert LSeyfert 2 ratio than is representative 
in the local Universe. 

In comparison to the active SMBH mass function con- 
taining the optically identified NL-AGNs of H04 (dotted 
line), we also find a significantly larger space-density of 
SMBHs. We find that the space-density of active SMBHs 
identified in the mid-IR is roughly constant in the mass 
region Mbh ~ (0.9-90) x 10 6 M with a value of $ w 
6.3 x 10" 4 Mpc" 3 logM"*. This space-density of AGNs is a 
factor of ~ 10 greater than that estimated by H04 over the 
same Mbh range. Since we find only two of the 17 AGNs 
in our sample are sufficiently luminous/unobscured to be 
detected in the SDSS survey, we determine that this is con- 
sistent with our results. We further suggest that the space 
density derived here may still be a lower-limit for the num- 
ber of NL-AGNs in the local Universe. A further examina- 
tion of the (distance-model adjusted) Palomar survey sug- 
gests that at least four further NL-AGNs are not included 
in our volume-limited survey. Of these, two (NGC 3486; 



11 The RBGS includes all IR-bright galaxies detected by IRAS 
with feoiim > 5.24 Jy at |6| > 5°. 



12 A. Goulding et al. 



I 



c 
o 



< 



00 



cn 
o 



u 

CL 



101 

cn 
o 



2 - 

3 7 
■4 L 
■5 
■6 
■7 h 



-8 r 
1.0 ~ 



0.5 - 



o.o - 



"rq 




■ D < 15 Mpc Active SMBH 

• Greene & Ho (2007) Broad-line SMBH 

Heckman et al. (2004) Narrow-line SMBH 

Marconi et al. (2004) Total SMBH 



H 1 1 h 



H 1 1 h 



H 

H 1 1 \ 



~19+{|% 



~18+) 2 % 

-*- 



+ 13°7 I- 



<14 + J. J % 



+ 



7 



8 



log(M BH ) [M ] 



Figure 5. (Upper panel) Comparison of volume-weighted space densities of active SMBHs in the local Universe, $ in units of number 
Mpc" 3 log Mgg. Mid-IR active SMBH function (filled squares; Goulding & Alexander 2009) is compared to the optically identified NL 
AGN function (dotted curve) of H04, the BL AGN function (filled circle) of Greene & Ho (2007), and total SMBH function (active+inactive 
galaxies; solid curve; Marconi et al. 2004). Sample selection bias is analysed using a robust Monte-Carlo simulation (shaded region; see 
Appendix). (Lower panel) Ratio of mid-IR active SMBHs to the total local SMBH mass function. The total SMBH mass function is 
extrapolated by 0.3 dex to Mbh < 10 6 Mq. For comparison the volume-weighted AGN fraction of H04 is also shown (dashed line). We 
estimate a mean volume-weighted local AGN fraction of 25+^ percent over the range M BH ~ (0.5-500) x 10 6 Mq. 



NGC 4565) lack high-resolution Spitzer-IRS spectroscopy 
of the central region (as noted in Table 2 of GA09), one 
(NGC 3031; L ra « 2.8 x 10 9 L ) lies fractionally below our 
luminosity limit for this survey, and NGC 4258 is not in- 
cluded in the RBGS due to its extremely large angular size. 

Whilst the main focus of this paper is to compare 
SMBH statistics derived from mid-IR and optical detec- 
tion techniques, it is prudent to note that the majority 
of AGN space densities calculated at higher redshifts are 
typically derived from sources detected in wide-field X-ray 
surveys. Hence, we now establish whether our mid-IR ac- 



tive space density may be missing a significant fraction 
of X-ray detected AGNs. Recently, a comparison between 
high-resolution Spit zer-IRS spect r oscop y and X-ray detected 
AGNs, was made bv lDudik etail (|2009l ) for a large sample of 
optically classified LINERS. Dudik et al. (2009) reported in- 
consistencies between [NeV] non-detections and the presence 
of hard X-ray nuclear emission in a subset of their sample. 
However, they conclude that the limited sensitivity of their 
mid-IR observations may be driving their observed result. 
Using the relations between high-ionisation mid-IR emission 
and hard X-ray luminosity (M08; GA09; Goulding et al. in 
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preparation), we suggest that a detected [NeV] luminosity 
of L[NeV] ~ 10 38 ergs _1 in an AGN (i.e., the limiting lumi- 
nosity in Dudik et al. 2009) would be equivalent to a hard 
X-ray luminosity of Lx,2-iokeV ~ 5 x 10 40 ergs~ 1 ; indeed, 
almost all of the X-ray detected AGNs which lack signifi- 
cant [NeV] emission in Dudik et al. (2009) are below this 
threshold. We thus conclude that with sufficiently sensitive 
high-resolution mid-IR spectroscopy, there is currently no 
conclusive evidence to suggest that X-ray detected Type-2 
AGNs lack significant [NeV] emission in the mid-IR. Hence, 
with the exception of some Seyfert 1 galaxies (as discussed 
previously) , it is unlikely that our derived mid-IR space den- 
sity lacks significant numbers of X-ray detected AGNs with 
[NeV] emission below our sensitivity limit. 

Given our comparatively small volume to that consid- 
ered by using the SDSS (e.g., Greene & Ho 2007; H04), we 
further validate our derived space density of active SMBHs 
by robustly testing our results to find if: 1) our sample is 
over-dense, and thus strongly subject to cosmic variance; 
or 2) given the modest errors associated with our Mbh es- 
timations, the derived space-density is strongly subject to 
scattering of objects in our defined binning structure. We 
discuss these analyses in the Appendix. Briefly, we find that 
our sample is broadly representative of galaxies to z ~ 0.3, 
and further show that even in our most pessimistic case, we 
find an increase in our derived space-density of at least a 
factor of ~ 2 (maximum increase by a factor of ~ 11) at 
M BH ~ 3 x 10 6 M Q over the optical NL space-density of 
H04. 



5.4 The volume-weighted local AGN fraction 

The ratio of the space densities of the active SMBH to total 
SMBH mass function (i.e., the volume- weighted local active 
SMBH fraction) is shown in the lower panel of Fig. [5] We 
calculate an overall active SMBH fraction of w 25^14 per- 
cent for SMBHs of Mbh ~ (0.5-500) x 10 6 M© down to our 
[NeV] completeness limit (LptoV] ^ 10 38 ergs _1 ). We find 
that this fraction is consistent with being constant through- 
out this Mbh range. However, given our detection sensitiv- 
ity limit, we are unable to probe lower Eddington ratios 
for AGNs hosting smaller SMBHs. Instead, we consider the 
effect of the AGN fraction for a fixed value of Eddington 
ratio (e.g., r\ > 10~ 3 ; i.e., thin-disc accretion systems) and 
find tentative evidence that the AGN fraction (w 16^g and 
~ 8^5° percent) may increase with decreasing SMBH mass 
(M BH « 10 6 -10 7 M Q and Mbh « 10 7 -10 8 M bins, respec- 
tively). 

For the lowest-mass SMBHs (Mbh ~ (5-30) x 10 5 M ), 
we estimate an overall non-negligible volume- weighted AGN 
fraction of 18ly 2 percent, potentially showing that a consid- 
erable proportion of small bulge (and pseudo-bulge) galaxies 
(i.e., late-type spiral galaxies; Sc-Sd) host AGN activity. It 
has been previously suggested by H04 and Greene & Ho 
(2007) that the AGN fraction may peak at Mbh ~ (0.7- 
2) x 10 7 M . However, with the inclusion of the additional 
low-mass optically unidentified AGNs (see Fig. [3J , we find 
that the AGN fractions are consistent with remaining con- 
stant or even increasing for Mbh < 10 7 M . 

As noted in §5.3, the space-density of AGNs, and hence 
the local AGN fraction derived in this work, may only be 
a lower-limit given the nature of our volume-limited survey 



which by definition does not include IR-faint systems. To im- 
prove upon these current source statistics, a larger sample 
of late-type spiral galaxies would be required to investigate 
our findings further. With the greater sensitivity and resolv- 
ing power proposed for the next generation of space-based 
mid-IR spectrographs, for example the Space Infrared Tele- 
scope for Cosmology and Astrophysics (SPICA l 12 l and the 
mid-IR instrument (MIRI) on-board the James Webb Space 
Telescope (JWST)^ surveys such as these can be continued 
and extended to study more distant (i.e., greater volumes) 
and heavily obscured AGNs. 



6 CONCLUSIONS 

We have presented the mean growth times and volume- 
weighted space density of active SMBHs in the local Uni- 
verse. Our sample of 17 AGNs was derived from a sensi- 
tive volume-limited mid-IR spectral survey of all IR bright 
galaxies to D < 15 Mpc carried out using the NASA Spitzer 
Space Telescope (see GA09 for further details on the sam- 
ple selection). The most accurate SMBH masses available 
for the objects are compiled from a variety of sources. For 
the three AGNs without published Mbh estimates, we use 
a bulge/disc decomposition method to determine the bulge 
luminosity and hence a SMBH mass (see section 3). Our 
main findings are the following: 

(1) Using combined mid-IR emission line and high-quality 
hard X-ray constraints, we have derived accurate measure- 
ments of the intrinsic luminosities of our sample of AGNs 
(see section 4). In conjunction with the well-established 
SMBH measurements from previous studies and our own es- 
timates from our aforementioned bulge/disc decomposition 
method, we have assessed the relative mass accretion rates 
of our sample. Due to our high sensitivity and the ability to 
probe low SMBH masses, we find that significant mass accre- 
tion (r) > 10~ 3 ) occurs onto SMBHs with Mbh « 1O 6 M , 
the majority of which would not be detected in even the 
most sensitive optical surveys. See section 5.1. 

(2) Using our derived relative mass accretion rates for 
the sample and assuming a typical accretion efficiency of e w 
0.1, we assessed the characteristic mean mass doubling times 
(fcjvf) for AGNs in the very nearby Universe. For AGNs 
hosting SMBHs with Mbh ~ (0.5-50) x 10 7 M we find 
consistent growth times (t%M ~ 47-198 Gyrs) with those of 
the NL-AGNs identified in the SDSS (H04). However, we 
also find that SMBHs with Mbh < 5 x 10 6 M (i.e., below 
the spectral resolution limit of the SDSS), are amongst the 
most rapidly growing SMBHs in the local Universe, with 
present-day growth times consistent with (and possibly less 
than) the current age of the Universe {p2M ~ 6-3 Gyrs). See 
section 5.2. 

(3) To assess the incidence of this population of low mass, 
rapidly growing SMBHs, we constructed a local space den- 
sity function of active SMBHs. We find that active SMBHs 
may be at least a factor ~ 2 more common than pre- 
viously identified in NL-AGN surveys using SDSS data. 
Furthermore, we estimate a non-negligible space density 



12 See 

13 See 



http: / /www. ir . isas . j axa.jp/ SPICA/ 



http://www.roe.ac.uk/uktac /consortium /miri/ 
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for low mass SMBHs (M BH ~ 1O 6 M ) of $ « 6 x 
10 -4 Mpc -3 log Mgg which is consistent with the space den- 
sity of more massive active SMBHs (Mbh ~ 1O 7 M0; i.e., 
those previously determined to be the most rapidly accret- 
ing population of SMBHs). See section 5.3 

(4) Using a local total SMBH mass function (Marconi 
et al. 2004), we estimate a mean volume- weighted local AGN 
fraction of ~ 23 jl^ percent, which remains relatively con- 
stant in the mass range M BH » (1-10) x 1O 6 M . However, 
when only considering the SMBHs with r\ > 10 - (i.e., radia- 
tively efficient accretion systems), we find tentative evidence 
for an increasing AGN fraction (» 16^g and « 8^5° percent) 
with decreasing SMBH mass (for M B h ~ 10 6 -10 7 M Q and 
Mbh ~ 10 7 -10 8 M , respectively). See section 5.4 
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APPENDIX A: VALIDATION OF DERIVED 
SPACE-DENSITY OF ACTIVE SMBHS 

As the space-density of active SMBHs derived in this work 
is significantly larger than that found in previous studies, in 
this section we attempt to validate our results by discussing 
possible limitations and additional sources of error which 
may exist in these analyses: 1) whilst the sensitivity of the 
data used in this survey is high, the volume considered is 
relatively small compared to that of the SDSS, thus our re- 
sults may be subject to cosmic variance; 2) given the modest 
errors associated with the Mbh estimates, the adopted Mbh 
binning structure is likely to be subjective and thus degen- 
erate towards objects scattering between the defined bins. 
Under these assumptions, in §A1 we investigate whether the 
sample is indeed representative of the local Universe and in 
§ A2 we discuss the construction of a Monte Carlo simulation 
to assess the effect of our adopted Mbh binning. 



Al Is our sample representative of the local 
Universe? 

Given the large incidence of AGNs within our sample it is 
possible that the volume considered in our sample is over- 
dense compared to other regions in the local Universe. The 
construction of the original sample of 64 bolometrically lu- 
minous galaxies in GA09 was designed to be complete down 
to the flux-limit of the Revised Bright Galaxy Sample. This 
imposed a distance constraint of D < 15 Mpc (see Fig. 1 
of GA09), and hence did not include the Virgo cluster at 
D ~ 16 Mpc, and thus our sample does not incorporate 
known local over-densities; however, this volume may not 
be representative of the Universe at large. 

To robustly test our considered volume (V ~ 1.3 x 

10 4 Mpc 3 ), we constructed a total SMBH space density 
function for all galaxies to D < 15 Mpc and compared 
this to the local total (active+inactive) SMBH mass func- 
tion of Marconi et al. (2004) derived from the luminosity 
function of local galaxies. Given the large co-moving vol- 
ume (V c w 1000 Gpc 3 ) considered in Marconi et al. (2004), 
their derived SMBH mass function is unlikely to suffer from 
significant cosmic variance. 

Our total SMBH mass function was formulated using 
all galaxies identified in the NASA/IPAC Extra-galactic 
Database (NED) to D < 15 Mpc with a total if-band lu- 
minosity of iK.gai > 1-5 x 1O 9 L0. The luminosity thresh- 
old is designed to include all galaxies which could poten- 
tially host a SMBH with Mbh > 10 6 M using the M B h- 
^K.bui relation. This conservative lower limit assumes that 
^K,gai = ^K.bui (i.e., that all galaxies in the sample have an 
early- type galaxy classification). In reality, the majority of 
the sources identified in NED to D < 15 Mpc are late-type 
galaxies (i.e., I/K, g ai 3> -£/K,t>ui), and therefore the Mbh limit 
is likely to be Mbh < 10 6 M Q . To D < 15 Mpc, we identify 

105 galaxies which potentially host a SMBH with M B h <5 

10 6 Mq. To estimate M B h in each of these galaxies, we re- 
late the associated Hubble-typ e from the Third Reference 
Catalogue of Bright Galaxies (Ide Vaucouleurs et al.l ll99lT ) 
to a mean bulge/disc ratio (e.g.. lBenson et alj|2007i ) and es- 
tablish an individual bulge luminosity based on galaxy-type 
and LK.bui- We convert the estimated bulge luminosity to 
M B h using the MBH-I/K,bui relation and construct a total 



SMBH mass function. Whilst this rather crude estimation 
carries large associated errors, we still find very good agree- 
ment (a mean variance of 0.1 dex) with the SMBH mass 
function of Marconi et al. (2004) throughout the mass range 
M BH ~ 10 6 -10 9 M . Thus, to first-order, our volume- limited 
sample does not appear to be over-dense and/or subject to 
strong cosmic variance, and hence, is broadly representative 
of the typical field-galaxy population in the local Universe. 

We further validate this conclusion by estimating the 
active SMBH mass function for those AGNs which are suf- 
ficiently optically bright to be included in H04. Under this 
assumption, the derived space densities of active SMBHs 
for this sample and that of H04 should be comparable if 
our sample is indeed representative of a field-galaxy popu- 
lation. We find that using the H04 detection limit inferred 
from Fig. [3] our new estimated space density is consistent 
with H04 ($ » 10 -4 ) throughout the mass range considered 
here (Mbh « (0.05-30) x 10 8 M Q ), although there are still 
considerable uncertainties given the small number statistics 
inherent with our sample. 



A2 Monte Carlo Analysis of Mbh binning 

We have established that our sample and the estimate of the 
SMBH mass density (both from active and inactive galax- 
ies) do not appear to be subject to over-densities caused by 
cosmic variance. We therefore now investigate the effect of 
small-number statistics on our results which are inherent in 
relatively small samples such as this. 

We show in §3 that accurate estimates of Mbh for our 
sample of AGNs from a homogeneous method is difficult 
to achieve. Thus, our adopted values of Mbh are derived 
from a variety of methodologies. In Fig. [S] the AGNs in the 
sample are placed into equal bins of Mbh with width 0.5 dex. 
However, this binning process does not allow for the error 
inherent to each individual Mbh measurement. Hence, some 
objects may scatter out of one defined bin and into another. 
The bin width we employ in our relatively modest sample (64 
objects) may therefore be subjective and requires testing. 
Here we use a Monte Carlo analysis to assess the effect of the 
scattering of AGNs into different Mbh bins on the derived 
space density of active SMBHs in the local Universe ($). 

Our Monte Carlo analysis calculates $ by selecting 
a random set of SMBH masses from gaussian probability 
distributions constructed using our adopted Mbh masses 
and their associated 1-a errors given in column 12 of Ta- 
ble 1. The binning structure was designed to incorporate 
at least one object in each bin from 250,000 realisations 
of the simulation across the considered mass range of the 
sample (Mbh ~ (0.5-500) x 1O 6 M ). We impose the upper 
mass limit as we have only one high-mass AGN (NGC 5128; 
Mbh ~ 2.4 x 10 s M©) in the sample. Thus, we determine 
that given our sample distribution the maximum number of 
equal width bins to be nine (i.e., equal bin sizes of 0.33 dex 
with i; 1 object). Hence, we use our simulation to conser- 
vatively assess the maximum error on our calculated space 
density of active SMBHs. 

The error shown in Fig. [5] (shaded region) is the stan- 
dard deviation of the 250,000 simulations in our volume- 
limited sample combined in quadrature with the Poisson er- 
ror determined from the number counts in our real sample. 
We found 250,000 realisations to be sufficient, since at this 



level the maximum variation in the standard deviation over 
multiple runs was less than 1CP 6 . 

We find that the mean spread in the derived value of 
$ from our volume-limited sample is w 0.85 dex, and as 
predicted appears to be subject to some scattering of Mbh. 
However, using our conservative error analyses, we show that 
the space density of active SMBHs found from our sample 
is consistently greater than that found for optical narrow- 
line AGNs (H04) with Mbh < 10 8 M . Thus, to first-order, 
we find a significant increase (of a minimum factor of ~ 2 
and a maximum of « 80) in the space density of active 
SMBHs in the local Universe in the mass range Mbh ~ (0.5- 
100) x 10 6 M . For Mbh > 10 s M our errors increase sig- 
nificantly due to very limited source statistics (see §4.1). We 
thus conclude that while this survey is subject to the scat- 
tering of objects through the Mbh bins, even in our most 
pessimistic case, we still find a significant increase in the 
space density of active SMBHs in the local Universe com- 
pared to that found in large-scale optical NL AGN studies. 
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